Solar Physics 
DOI: 10.1007/' 



Magnetic Field Structures in a Facular Region 
Observed by THEMIS and Hmode 

Y. Guo^^ • B. Schmieder^ • V. Bommier^ • 
S. Gosain^ 

© Springer •••• 

Abstract The main objective of this paper is to build and compare vector 
magnetic maps obtained by two spectral polarimeters, i.e. THEMIS/MTR and 
Hinode SOT/SP, using two inversion codes (UNNOFIT and MELANIE) based 
on the Milne-Eddington solar atmosphere model. To this end, we used observa- 
tions of a facular region within active region NOAA 10996 on 23 May 2008, and 
found consistent results concerning the field strength, azimuth and inclination 
distributions. Because SOT/SP is free from the seeing effect and has better 
spatial resolution, we were able to resolve small magnetic polarities with sizes of 
1" to 2", and we could detect strong horizontal magnetic fields, which converge or 
diverge in negative or positive facular polarities. These findings support models 
which suggest the existence of small vertical flux tube bundles in faculae. A 
new method is proposed to get the relative formation heights of the multi-lines 
observed by MTR assuming the validity of a flux tube model for the faculae. We 
found that the Fe i 6302.5 A line forms at a greater atmospheric height than the 
Fe I 5250.2 A line. 

Keywords: Active Regions, Structure; Magnetic fields. Photosphere; Polariza- 
tion, Optical 



1. Introduction 

Faculae are bright plages observed in Ha or photosphcric lines. They are consid- 
ered as magnetized regions constituted of a bundle of thin vertical fiux tubes with 
a magnetic field strength of thousands of gauss and a size of tens to hundreds of 
kilometers ( |Zwaan, 1987p . Similar values for field strength and size of the flux 
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tubes have been found by analyzing the ratio of Stokes V signals in the Fe i 5250.2 
A, Fe I 5247.1 A, and Fe i 5232.9 A lines, for a spatially unresolved photospheric 
network ( [Stenflo, 1973[ ) . The flux tubes are located at the boundaries of granules 
and bundled by the convection flow. Since the pressure of the solar atmosphere 
decreases with height, the tubes expand and the magnetic field strengths decrease 
to compensate the pressure decrease. The network field expands rapidly with 
height and forms a canopy above the internetwork region. Similarly, the faculae 
would be bundles of flux tubes expanding with height. 

In order to derive the three-dimensional (3D) magnetic structure of these 
features, we require multi-line spectral polarimetry. On the other hand, to re- 
solve the fine structure, we require high spatial resolution. The French-Italian 
ground-based solar telescope THEMIS (Telescope Heliographique pour I'Etude 

du Magnetisme et des Instabilites Solaires) on the Canary Islands ( [Lopez Ariste, Rayrole, and Semel, 2000H 
Bommier et al., 2007[ ) and the Japanese space borne satellite Hinode ( Kosugi et al., 20"07H 



provide this ability. THEMIS is designed to be free from the instrumental polar- 
ization, since the polarization analyzer is located at the primary focus. It is also 
characterized by the multi-line spectroscopy capability in the MTR (Multi-Raies) 
grid mode that ensures the best co-spatiality. The beam exchange technique of 
MTR increases the polarimetric accuracy as the flat field errors are minimized. 
Solar Optical Telescope (SOT) aboard the Hinode spacecraft is not affected 
by the Earth's atmospheric seeing effect. It observes diffraction limited im- 
ages with 0.2" - 0.3" resolution by its 50 cm aperture optical telescope. The 
Spectral Polarimeter (SP) behind SOT observes the full Stokes profiles of two 
magnetically sensitive lines Fe i 6301.5 A and Fe i 6302.5 A (jlchimoto et al.\ 
IShimizu et dX [20051 ISuematsu etal\ [MM ITsuneta et al\ [MI5|) . With 
its high performance, breakthroughs in different topics related to photospheric 
fields have been obtained by SOT: the detection of horizontal fields over granules 
and the hidden turbulent magnetic flux ( Lites et al., 2007 ), the birth of small 
flux tubes with kilo- gauss field strength ( [Nagata et al., 2008[ ), and the discovery 
of transient horizontal magnetic fields in quiet sun Pshikawa et al., 2008[ ). 

The aim of this paper is: (i) to demonstrate that the two spectral polarimeters 
THEMIS/MTR and Hinode SOT/SP, and the two Stokes profile inversion codes 
give consistent results, and (ii) to study the magnetic structure of faculae. We 
observe a facular region, build and compare vector magnetic fields using the 
full Stokes profiles of multi spectral lines observed by two spectral polarimeters 
(MTR and SOT/SP) and fitted by two inversion codes. We find concentrations 
of magnetic fiux with converging and diverging transverse field vectors in high 
angular resolution magnetograms observed by SOT/SP. Finally, we use a flux 
tube model to get a qualitative view on the formation heights of two Fe i lines. 
The multi-line capability of THEMIS/MTR and the high angular resolution of 
SOT/SP complement each other in this study. The description of the observa- 
tions and the methods of data analysis are given in Section [5] The results are 
presented in Section [31 We discuss the magnetic field gradient according to the 
fiux tube model in Section [31 The conclusions arc drawn in Section [SI 
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2. Observations and Data Analysis 
2.1. Instruments and Observations 

The faculae in the active region NOAA 10996 (N10° W32°) were simuhaneously 
observed by THEMIS and Hinode. THEMIS/MTR scanned the region from 
16:44-17:38 UT and Hmode SOT/SP from 14:30-15:01 UT, on 23 May 2008. 
The hehocentric angle for the observed region was « 34°. We plot the full field 
of view of SOT/SP and THEMIS/MTR in Figure ffl Both instruments scan the 
solar surface from east to west. THEMIS/MTR recorded the two beams on the 
same CCD camera using a grid in order to avoid problems of co-spatiality. The 
grid had empty space (15.5") between the bars, along the slit (« 120"). The grid 
bars were covered with three 8" steps along the slit (which was north-south). 
The next step was to move the telescope towards the west with a step size of 
0.8". The scan of the region (110" x 104") was obtained within about one hour. 
The characteristics of THEMIS/MTR are given in Table [H 

The three photospheric lines observed with MTR were selected because they 
are all normal Zeeman triplet lines (formed by transition between energy levels 
with J = and J = 1). Moreover, their Lande g factors are rather large, 
3.0, 2.5 and 2.0 for Fe i 5250.2 A, Fe i 6302.5 A, and Ca i 6102.7 A, respec- 
tively. For both reasons, these lines in the visible spectrum have good magnetic 
sensitivity. They are also well adapted to the UNNOFIT inversion, which can 
treat only normal Zeeman triplet lines. For the other lines we have UNN0FIT2 



(Bommier et ai, 2007), but generally their Zeeman sensitivity is lower and the 
inversion fails to converge in pixels with poor signal-to-noise ratio. 

The fast map mode of SOT/SP provides us with a relatively high cadence of 
3.8 s per position of the slit, and a spatial resolution of « 0.3" per pixel after 
binning on flight. The width of the slit was 0.16", and the length was reduced 
to 164"/2 for this study. The scan of the region (145" x 82") was obtained 
within 30 minutes. We selected a central part of the active region, where data 
were observed with good quality both by SOT/SP and MTR (solid rectangle in 
Figure [T]), to perform a detailed analysis. 

2.2. Inversion Procedures and Bisector Method 

How can we extract the vector magnetic field from the full Stokes profiles? It 
needs two steps. The forward problem answers how the line forms for given mag- 
netic fields in the photosphere of the Sun, by solving the radiative transfer equa- 
tions for polarized radiation ( |Unno, 1956[|Rachkovsky, 1962[|Jefferies, Lites, and Skumanich, 1989[ ).| 
The inversion problem extracts the magnetic field and other physical parameters 
by fitting the observed line profiles with the theoretical line profiles obtained in 
the first step ( Landolfi and Landi Degl'Innocenti, 1982[ Skumanich and Lites, 1987] ) .[ 



All the raw spectra of THEMIS are calibrated by spectral destretching, dark 
current subtraction and flat field correction ( [Bommier and Rayrole, 2002[ [Bommier and Molodij, 2002p | 
before polarimetric analysis and inversion. In Table [1] we list the observational 
characteristics of five of the major lines that we have used within the wavelength 
range of MTR, and of the Fe i 6302.5 A line of SOT/SP. The I,Q,U, and V 
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profiles of SOT/SP are obtained by adding and subtracting tlie raw spectra 
after calibration with the standard data analysis packages for SOT/SP in Solar 
Software (SSW) Pchimoto et al., 2008[ ). 



We use the UNNOFIT inversion code ( Landolfi and Landi Degrinnocenti, 1982H 



[Bommier et al., 2007[ ) to fit the line profiles of Ca i 6102.7 A, Fe i 6302.5 A, 
and Fe i 5250.2 A observed by THEMIS/MTR. The data obtained by SOT/SP 
are fitted by two inversion codes, i.e. by UNNOFIT and by the code devel- 
oped by Hinode team based on MELANIE (Milne-Eddington Line Analysis 
using an Inversion Engine; ISocas-Navarrol [M?T|) . UNNOFIT and MELANIE 
use the Levenberg-Marquardt algorithm to achieve the least fitting of the 
observed profiles by the model profiles. The model profiles are given by the Unno- 
Rachkovsky solution based on the Milne-Eddington approximation for the solar 
atmosphere. Including the magneto-optical and damping effects, all parameters 
needed to synthesize the theoretical Stokes profiles are listed in Table [51 The 
damping parameter 7 and damping constant a arc related by a = 7/(47rAAD)- 
The parameter of the source function b is defined hy b — ^Bi/Bq, where fi is 
the cosine of the angle between the normal to the solar surface and line-of-sight 
(LOS), and B{t) = Bq + Bit, where r is the optical depth. IBommier et al.\ 
p007|) introduced the fiUing factor in UNNOFIT as follows: 



(1) 



Um, and Vm are the Stokes parameters in magnetic fields, /q is the 
intensity of the quiet sun. They have the same physical conditions except for the 
presence or absence of magnetic fields. So the profiles of /q in UNNOFIT are cal- 
culated using the radiative transfer equations with the same parameters as used 
for /m, while ignoring the magnetic fields. Note that /q in MELANIE is given by 
the average intensity profiles within the whole field of view excluding active re- 




gions, and it is usually recognized as the stray light (Skumanich and Lites, 1987), 
whose fraction / is relative to the filling factor a as a = 1 — /. All the magnetic 
field maps that are plotted in the figures of this paper have been multiplied 
by the filling factors. They are aB for UNNOFIT inversions and (1 — f)B for 
MELANIE inversions. These local averaged magnetic fields are the only deter- 
mined quantities, since the intrinsic field strengths and filling factors cannot be 
determined separately as shown bv IBommier et aLI(|2007p . From now on, we will 
use B instead of aB to represent the local average field strength for concision. 
But we should keep in mind that all field strengths and field components, such 
as the LOS magnetic fields, have been multiplied by the filling factors. 

The magnetic fields of the Na Di line profiles observed by THEMIS/MTR 
are calculated by the formula ( [Semel, 1967[ ) as follows: 

SX = 4.67 X 10-"5^ff a^Slos , (2) 

where g^s = 1.33 for the Na Di line, and the units of wavelength and mag- 
netic field arc A and gauss, respectively. The formula shown in Equation ^ is 
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different from the weak field approximation, where the LOS magnetic field is 
calculated hy V oc BLOsfj- Equation ^ remains valid even for strong fields 
( [Bommier, Rayrole, and Eff-Darwich, 2005] ). We determined the Zeeman shift 
5X by a bisector method using the middle point of a chord with a length of 
2AA, which connects the equal intensity points of a line profile. The wavelength 
difference of the middle points of the I + V and I — V profiles is the Zeeman shift 
6X dBerlicki, Mein, and Schmieder, 2006t , and A A represents the wavelength dif- 
ference from the line core approximately. We obtain the magnetic fields for AA 
= 0.08 A, 0.16 A, 0.24 A, and 0.32 A in the Na Di line. 



2.3. Co-alignment 



In order to compare the magnetic field characteristics, we have to co-align 
the images in different wavelength bands and by different instruments. The 
observations in different wavelength bands of MTR have offsets between each 
other. The pixel sizes for these wavelength bands also differ slightly, which 
requires an interpolation to the same pixel resolution before co-alignment. Also, 
the instruments observe only partial regions of the solar disk; hence pointing 
errors cannot be avoided. We use the magnetic map of the Michclson Dopplcr 
Imager (MDI; IScherrer et ail I1995P as the basic reference frame, since the full- 
disk observation ensures the accuracy of the MDI coordinates by comparing the 
position of solar limbs ( |Guo et ai, 2009] ). Firstly, we make a coarse alignment of 
the MTR and SOT/SP magnetic maps referred to the MDI magnetic map that 
has the closest observation time. The so-called feature identification method is 
employed, where several points sharing common features on the two images are 
selected consecutively to determine the shifts between them. Next, the magnetic 
maps of MTR are interpolated to a common spatial resolution, which is that of 
Fe I 6302.5 A here. The magnetic map of SOT/SP is differentially rotated to the 
middle observation time of the MTR scanning and interpolated. Finally, the fine 
alignment is done by correlating each image with the Fe i 6302.5 A magnetic 
map, and finding out the position where two images have the largest correlation 
coefficient. 

We present in Figure [2] some of the co-aligned maps of the LOS magnetic 
fields in the selected field of view, shown by the solid rectangle in Figure [T] For 
MTR, magnetic fields have been calculated for the Ca i 6102.7 A, Fe i 6302.5 
A, Fe I 5250.2 A, and Na Di 5896 A lines, and for the SOT/SP Fe i 6302.5 A 
line. For the Na Di line we have calculated four magnetic maps at AA = 0.08 
A, 0.16 A, 0.24 A, and 0.32 A. Only one Na Di magnetic map for AA = 0.08 
A is shown Figure [5] as an example. We draw two cuts in each magnetic field 
map in order to perform a detailed comparison of the field distributions in the 
subsequent sections. 



2.4. 180° Ambiguity Removal 



The transverse component of the magnetic field vector has an intrinsic 180° 
ambiguity. We have to remove this ambiguity before we can transform the 
magnetic field to the local heliographic coordinates for a better understanding 
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of the magnetic field structure. IMetcalf et aL \ (|2006p reviewed different existing 
methods that aim to remove the f80° ambiguity. A new method based on the 
fields measured at two different heights and the divergence-free nature of the 
magnetic field has been recently proposed by [Crouch, Barnes, and Leka| (|2009p 
and is currently under development. 

fn this work we adopt the method of non-potential magnetic field calcula- 
tion, developed by |Georgoulis| ((2005|) and IMetcalf et aL\ (|2006| . This method 
iteratively compares the observed magnetic field vectors with the model field 
B = Bp -|- Be in the heliographic coordinates system, where the potential com- 
ponent Bp is calculated by potential field extrapolation and the non-potential 
component Be is calculated by the direct and inverse Fourier transforms of 
the vertical electric current density ( |Chae, 2001[ |Georgoulis, 2005] ) . We use 
the improved version of the non-potential magnetic field calculation method 
( [Metcalf et ai, 2006[ ), where the vertical electric current density Jz is updated 
in every iterative step based on the intermediate ambiguity solution B. is 
calculated by means of Ampere's law. The field of view for the potential and 
non-potential magnetic field is padded with zeros to allow for a more accurate 
implementation of the Fourier transform. 



2.5. MTR Ha Observation 



We plot the intensity map of the facular region at the line center of Ha 6563 A in 
the bottom panel of Figured The Ha map is not precisely aligned with the mag- 
netic field map, but the alignment is sufficient for a qualitative comparison. The 
comparison shows that dark filaments are located between positive and negative 
magnetic polarities. The Ha intensity increases in the regions of magnetic field 
concentration. This can also be found in the Hinode SOT/BFl (Broadband Filter 
Imager) Ca ii H line intensity map. It is still not fully understood why taeniae 
arc brighter than a normal quiet sun region. One possible explanation is that the 
optical thickness is smaller in magnetic tubes than in the external photosphere. 
This effect allows us to see deeper and hotter layers ( [Strous, I994P . IKeller et al.\ 
(|2004p have performed a 3D simulation of non-grey radiative magneto-convection 
to explain the origin of this phenomenon. [Hasan and van Ballegooijen| (I2008P 
suggest that the bright features are heated by the dissipation of short period 
magneto-acoustic waves in magnetic flux tubes. 



3. Results 

3.1. Magnetic Field Maps and Histograms 

Magnetic maps of MTR Ca i 6102.7 A, Fe i 6302.5 A, Fe i 5250.2 A and SOT/SP 
Fe I 6302.5 A are shown in Figure [S] All magnetic field parameters have been 
transformed into a local solar coordinate system. The arrows are plotted for 
the transverse fields with strength larger than 50 gauss, which is the limit of 
the transverse field accuracy ( [Ichimoto et al., 2008p . The results obtained from 
MTR show that in the regions of strong concentrations, the magnetic fields are 
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almost vertical. The horizontal components are typically less than 50 gauss; 
therefore, only a few arrows appear in these regions. In the central part of the 
facular region the arrows could correspond to the footpoints of Ha fibrils that 
appear in the bottom panel of Figure [2] between the two bright regions. 

The fields obtained by SOT/SP show more detailed features. The horizontal 
vectors at the inner edges of the strong field concentration regions also indicate 
that this region could correspond to footpoints of the Ha fibrils. In large concen- 
trations we can detect many small cells with horizontal fields pointing towards 
and away from the center of negative and positive polarities, respectively; sec the 
circles in the bottom panel of Figure [31 This indicates the existence of expanding 
flux bundles rooted in these cells, supporting the canopy model. 

The histograms of the field strength, azimuth, and inclination angles that 
have been transformed into the heliographic coordinates after the 180° ambiguity 
removal are shown in Figure 3] There are two peaks that appear in each of the 
histogram of the azimuth angle. The plots are in the range of —180° to 180°, 
since the ambiguity has been removed. The histograms of the inclination angles 
obtained from MTR Fe i 6302.5 A and Fe i 5250.2 A clearly show that the 
magnetic fields are mainly vertical in the local solar surface plane, and the incli- 
nation histograms of SOT/SP Fe i 6302.5 A show departures from the vertical 
direction. This is due to the expansion of the flux tube and the high spatial 
resolution of SOT/SP. MTR observes the average value indicating mostly the 
vertical component. 

The inclination angles measured by the THEMIS/MTR Fe i lines are more 
vertical than that measured by the Ca i 6102.7 A line. There are several possible 
explanations for this difference. First of all, the uncertainties of the angle mea- 
surements may account for it. IBommier et al.\ (|2007p found that azimuth and 
inclination angles can only be determined within the accuracy of 10° to 20°, due 
to the polarimetric accuracy of THEMIS/MTR. Secondly, the magnetic flelds 
obtained by different lines are formed at different altitudes and the Ca i line 
could be formed at greater atmospheric heights than the other lines, where the 
magnetic fleld is more inclined. Finally, the Ca i 6102.7 A line may not be formed 
in local thcrmodynamical equilibrium (LTE), whereas UNNOFIT is based on the 
Milne-Eddington model of the atmosphere which assumes LTE. In that case the 
magnetic flelds obtained from the inversion of the Ca i 6102.7 A line would have 
large uncertainties. 

3.2. Comparison of Magnetic Fields Observed by THEMIS/MTR and SOT/SP 

Hinode observes seeing-free images with better quality and provides flner struc- 
tures in this facular region than ground-based instruments. It makes sense to 
compare the data from Hinode SOT/SP and THEMIS/MTR. In order to make 
the comparison reliable, we treated the Stokes proflles of SOT/SP with both 
UNNOFIT and MELANIE. The comparison results are shown in Figure El The 
magnetic flelds obtained by UNNOFIT and MELANIE inversion codes are found 
coherent, especially for LOS components. Yet, the absolute values obtained by 
UNNOFIT inversion code are generally larger than that by MELANIE inversion 
code, which can be explained by the different free parameters considered in the 
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two inversion codes, especially the different methods of the computation of the 
filling factors. Only a few points violate the positive correlation relationship, and 
most of them appear in the weak magnetic field regions. 

The structures are clearly shown in the top row of Figure |6l where the LOS 
magnetic fields obtained by SOT/SP and MTR Fe i 6302.5 A are plotted along 
two selected cuts. The absolute value of SOT/SP LOS magnetic fields reaches 
about 1200 gauss, while the one of MTR is about 300 gauss. The sizes of the 
structures are about 1" to 2" with SOT/SP if we define them as the full- 
width at half-maximum of the absolute value of the field strength variation 



along a cut, which are comparable to previous direct measurement ( |Miiller, 1994 



Daras-Papamargaritis and Koutchmy, 1983). These values should correspond to 



the sizes of the bundles of flux tubes. The sizes determined by the observation 
of MTR in the same region are larger. The magnetic fields observed by SOT/SP 
are obtained by the MELANIE inversion. We have already shown in Figure 
[S] that MELANIE and UNNOFIT give consistent inversion results. However, 
considering that MELANIE has been optimized for SOT/SP data by the Hinode 
team, we still use the MELANIE inversion result for the SOT/SP data to be 
compared with the UNNOFIT inversion result of MTR data. The time difference 
between the two observations is about 2 hours. We have checked the evolution of 
this region by four sets of magnetograms taken on the same day, and we found 
no evidence for a decreasing trend in the field strengths. 

The difference of field strengths can be understood by the difference of pixel 
sizes (0.3" for SOT/SP and 0.8" for THEMIS/MTR),and by considering the 
seeing effects and the photometric accuracy. THEMIS observes on the ground, 
where the atmosphere of the Earth blurs, distorts the observed images, and 
decreases the spatial resolution. We adopt a convolution process to simulate the 
seeing effect on the Stokes profiles observed by SOT/SP. Each Stokes profile is 
reconstructed after a convolution of each wavelength image by a 3" x 3" two- 
dimensional (2D) boxcar function sampled at the SOT/SP pixel size. The photo- 
metric accuracy of MTR is « 3 x 10"'^ ( Lopez Ariste, Rayrole, and Semel, 2000 ) 
and it is w 10"^ for SOT/SP Pchimoto et al., 2008[ ). We add a random noise 
with a normal distribution in the field of view and with a standard deviation 
of 3/c X 10"'^ to each Stokes profile of SOT/SP at each wavelength, /c is the 
continuum intensity. Then the MELANIE inversion is applied. The distributions 
of the LOS magnetic fields on the selected cuts are shown in the bottom row 
of Figure \E[ Both the field strength and the size of the magnetic concentration 
feature of the degenerate SOT/SP LOS magnetic field map are comparable to 
what MTR observed. 



3.3. Formation Heights of Ca i and Fe i Lines 



Some simulated formation heights of Fe i lines are listed in Table [31 The heights 
are measured from T5000 = 1, and T5000 is the optical depth at the wavelength 
of 5000 A. All the calculations adopt the LTE assumption, but with different 
atmosphere models. [Bommier et aLI(|2007p employed the quiet sun photospheric 
reference model ( [Maltby et al., 1986[ ), extrapolated downwards beyond —70 km 
to —450 km below the T5000 = 1 level. Above —70 km, this model is very 
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similar to the quiet sun VAL-C model ( [Vernazza, Avrett, and Loeser, 198T |. 
|Bruls, Lites, and Murphy | p99ip adopted the VAL-C model, and ISheminoval 
(|1998p assumed the Holweger-Miiller model ( |Holweger and Miiller, 1974D . The 
formation heights show some discrepancies. The calculation of these heights 
depends on the choice of the physical parameters, such as the temperature, 
magnetic field, and microturbulent velocity. We would need the formation al- 
titudes of lines by using the response function to magnetic field perturbations 
(Vitas, private communication). It is not in the scope of this paper. However, 
we propose another way to derive the relative formation altitudes of the lines by 
using the assumption of the flux tube model in the faculae. 

In order to have a quantitative comparison of the LOS magnetic fields ob- 
tained by the different lines of THEMIS/MTR, we produce scatter plots, shown 
in Figure [71 The plots show that for each pixel the LOS magnetic fields obey the 
following relation, except for some weak field pixels that we will not consider, 
because they arc not in facular concentrations: 



|SlOS,6102.7| > I^LOS, 6302.51 > |-SloS, 5250.2 | • (3) 

Presently, we have measured the magnetic field strengths in different lines. In 
the following sections we will use the LOS magnetic field, but not the strength, 
for two reasons. First, the field in the region is mainly vertical and the region is 
not so far from the solar disk center. Secondly, only the LOS magnetic field can 
be computed for the Na Di line. In the flux tube model the tube expands and 
the magnetic field strength decreases with height. Equation ^ would indicate 
that the formation height of Fe i 5250.2 A is the highest one, and the Ca i 6102.7 
A is formed below the Fe i 6302.5 A. 

We can also test this assumption of the flux tube model by measuring the 
full-width at half-maximum of the absolute value of the field strength variation 
along a cut in the image. We analyzed the LOS magnetic fields along the selected 
cuts presented in Figure [5J They are plotted in the top row of Figure [51 The 
full-widths at half-maximum follow the relations: 



W^LOS,6102.7 > W^LOS,6302.5 > W^LOS,5250.2 • (4) 

In the flux tube model we expect that the structure with the smallest width 
would be formed at the lowest height. Equation ([4]) is in contradiction with what 
Equation ([3]) implies. We can use the measurements of the magnetic field along 
the Na Di line wings, where the relative formation heights of each wavelength 
have already become known, to explain this inconsistency. 



3.4. Magnetic Field Gradient in the Na Di Line Wings 



We have computed the fields at different wavelengths in the line wing of Na 
Di 5896 A, i.e. AX = 0.08 A, 0.16 A, 0.24 A, and 0.32 A. The LOS magnetic 
fields along the white and black cuts as shown in Figure [H are plotted at the 
bottom row of Figure [51 We measured the full- widths at half-maximum of the 



SOLA: ms.tex; 23 February 2010; 15:41; p. 9 



Guo et al. 



structures, which peak at 6 to 12 arcsec in all the plots of FigureEl and we found 
the following relation: 

W^LOS,5896±0.08 > W^LOS,5896±0.32 • (5) 

The magnetic features observed in the line core have a larger size than that in 
the line wing. This result is consistent with the flux tube model assumption. Let 
us analyze the magnetic field strength variation along the Na Di profile. 

From the bottom row in Figure [51 we got the following relation of the mag- 
netic fields obtained in the core and wing of Na Di (AA = 0.08 A and 0.32 A 
respectively): 

|-BlOS,5896±0.08| > |-SlOS,5896±0.32| , (6) 

The fields in the line core of Na Di are stronger than those in the line wing. We 
found a similar contradiction by using Ca i and Fe i lines. Such a behavior for 
the Na Di line has been reported by |Berlicki, Mein, and Schmieder (j2006p and 
IMein et aLI([M?7| . 



4. Discussion 

4.1. Magnetic Flux Tube Modelling 

IMein et al.\ (j2007p resolved this contradiction by a 2D magnetic flux tube model. 
They computed theoretical profiles of the Na Di line in an atmosphere where 
fiux tubes and quiet sun regions coexist side-by-side in a pixel. First, they showed 
that the increase of magnetic field strength with height can be explained by a 
pure geometrical effect at the edge of flux tubes. For one pixel with a size of 
2s, where s is the half- width of the spatial resolution, a smaller quiet sun region 
should be considered at higher levels than at lower altitudes {left panel of Figure 
E]). Secondly, it should be explained by the combination of small filling factors 
with slopes of the Stokes profiles at the flux tube center {right panel of Figure 
In fact, the measurements of Stokes parameters and magnetic field (/, V , 
and B) are obtained by THEMIS with a spatial resolution, 2s, much larger than 
the expected size of the fiux tube sections 2wA\. [Mein et an(|2007p studied this 
effect on the Stokes profiles by using flux tube models with different sizes. On 
the one hand, the filling factor of the lower formed line is smaller because the 
size of the flux tube is smaller. On the other hand, the Stokes / profile slope 
depends on the atmosphere model. In the line wing, the absolute value of the 
profile slope ^ is larger in the quiet sun than that in the flux tube center. It is 
reversed in the line core. Assuming that /q(AA) and /m(AA, 0) are the intensities 
in the quiet sun region and magnetic flux tube center as denoted in Figure [9j we 
have the following relationship for |AA| >0.24 A: 

a/q(AA) a/.n(AA,0) 

' aA ' ' d\ ' ' ^ ' 
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and for AA <0.24 A: 

— dx — 

The ratio 'r(AA, 0) of the profile slope between the quiet sun and flux tube region 
is larger in the line wing than in the line core, where 

IMein et al.\ (|2007p derived the following equation from a simpler but compre- 
hensive model, where the magnetic field in the fiux tube is constant at the same 
height: 

The half-widths of the tube and the spatial resolution of the instrument are 
denoted by waa and s, respectively, as shown in Figure [HI i?m(AA,0) and 
i3(AA, 0) are the LOS magnetic field in the flux tube center and the synthetic 
magnetic fleld convolved with the quiet sun region within the spatial resolution 
2s, respectively. B{A\, 0) corresponds to the observed magnetic field. 

i3(AA, 0) has a negative relationship with the ratio r(AA, 0), which is larger 
in the line wing, and a positive relationship with waa, which is smaller in the 
line wing. So i?(AA,0) will become smaller at lower altitudes after convolving 
the profiles of the flux tubes with the profiles originating from the quiet sun. 

The discussion above is for a single flux tube, whose size may be smaller 
than the spatial resolution of SOT/SP (0.3"). In this case, each of the peaks 
Al, A2, Bl, B2, and B3 measured by SOT/SP shown in Figure IH] is constituted 
by a bundle of smaller flux tubes. IMein et al.\ (|2007p showed that the sizes of 
the flux tube bundle still increase with height in a multi-flux tube simulation. 
Further, Equation (10) still holds, i.e. a combination of Stokes proflles in quiet 
and active regions may reverse the field strength variation with height. 

4.2. Relative Formation Heights of the Fe i, Ca i, and Na Di Lines 

The results of Equation in contradiction with the half- width relation, are 
now explained by the small filling factor of the magnetic field and the difference 
of slopes of the Stokes / profiles observed in quiet sun and in the flux tube. We 
believe that these effects will also affect the Ca i and Fe i lines. Thus, we suggest 
that the formation altitude of Ca i 6102.7 A is higher than that of Fe i 6302.5 
A, which forms higher than Fe i 5250.2 A. The formation height of Ca i 6102.7 
A is comparable to the one where the line wing of Na Di forms. Nevertheless, 
we should still keep in mind that the magnetic fields obtained by Ca i 6102.7 A 
were computed in LTE and the Milne-Eddington approximation, which may be 
not adequate for this line. The results for the two Fe i lines are more reliable; 
however, the formation height of a line varies along the line profile itself. In 
addition, the magnetic inversion makes use of the whole profile, which forms at 
different altitudes. Nevertheless, if Fe i 6302.5 A forms at greater heights than Fe 
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I 5250.2 A at each wavelength from the hne center, this relation will be conserved 
by the inversion. 



5. Conclusion 

In this paper we have compared the vector magnetic fields obtained from the 
THEMIS/MTR muhi-hne observations and Hinode SOT/SP observations for 
a facular region. Two inversion codes, i.e. UNNOFIT and MELANIE, were 
adopted to determine the field strength, azimuth, and inclination angles in this 
region from the full Stokes profiles. The 180° ambiguity was removed by the non- 



potential magnetic field calculation method ( Georgoulis, 2005 Metcalf et al., 20()6[ ) 
The magnetic fields in the faculae are found to be mainly vertical. We resolve 
small flux tube bundles in the faculae with converging and diverging horizontal 
components of magnetic flelds in negative and positive polarities, respectively. 
The maximum field strengths observed by SOT/SP (1200 gauss) are four times 
larger than those obtained from MTR. The sizes of the detectable facular details 
are found to be about 1" to 2" for SOT/SP, while they are larger when measured 
by MTR. Nevertheless, we conclude that MTR and SOT/SP give consistent 
results concerning the characteristics of the fields, particularly for the inclination 
and azimuth histograms. The differences between the field strengths are reduced 
to a factor of « 1.5 if the Earth's atmospheric seeing effects are taken into 
account. 

The assumption of the conservation of magnetic flux in the expanding flux 
tube model has been used to derive the formation heights of the lines observed 
with MTR. We also consider the opposite effect on the magnetic field strength 
gradient that is generated by poor spatial resolution compared to the size of the 
fiux tube dMein et al, 2007D . Fe i 6302.5 A forming at greater heights than Fe 
I 5250.2 A can be compared with some simulated results listed in Table [3] It 
is consistent with the calculation by IBommier et al.\ (|2007p obtained by solving 
the radiative transfer equation. 

The combination of the high spatial resolution of SOT/SP and the multi- 
line observation of THEMIS/MTR allowed us to have a 3D view on the flux 
tubes in faculae. To get a quantitative model, one would need to calculate more 
complete and reliable formation heights of the lines using response functions. 
The formation height of different magnetically sensitive lines is an important 
parameter in the new 3D method to resolve the 180° ambiguity, where the 
divergence-free nature of magnetic fields is used. Removing the 180° ambiguity 
is necessary to built vector magnetic fields, which are used as the boundary 
condition for nonlinear force-free field extrapolations to get the fields in the 
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Figure 1. The full fields of view of SOT/SP (dashed rectangle) and THEMIS/MTR (dotted 
rectangle). The solid rectangle denotes the region that we select to study the magnetic field in 
detail. The background is part of the full disk 96 min magnetogram observed by SOHO/MDI. 



Table 1. Instrumental characteristics for different lines of THEMIS/MTR 
and Hinode SOT/SP. The data were observed on 23 May 2008. 



Line 


Scan duration 


Dispersion 


Pixel size 




(UT) 


(A/pixel) 


(arcsec) 


MTR Ca I 6102.7 A 


16:44:03 - 17:38:55 


0.0130 


0.800, 0.236 


MTR Fe I 6302.5 A 


16:44:03 - 17:38:55 


0.0125 


0.800, 0.233 


MTR Fe I 5250.2 A 


16:44:03 - 17:38:55 


0.0111 


0.800, 0.230 


MTR Na Di 5895.9 A 


16:44:03 - 17:38:55 


0.0130 


0.800, 0.232 


MTR Ho 6562.8 A 


16:44:03 - 17:38:55 


0.0144 


0.800, 0.226 


SOT/SP Fe I 6302.5 A 


14:30:05 - 15:01:38 


0.0215 


0.297, 0.320 



P. Mein, M.K. Georgoulis, N. Vitas, and T. Torok very much for their help and discussions. 
Y. Quo was supported by the scholarship granted by China Scholarship Council (CSC) under 
file No. 2008619058. S. Gosain acknowledges CEFIPRA funding for his visit to Observatoire 
de Paris, Meudon, France under its project No. 3704-1. 

References 

Berlicki, A., Mein, P., Schmiedcr, B.: 2006, Astron. Astrophys. 445, 1127. 



250 



200 



150 



100 



SOLA: ms.tex; 23 February 2010; 15:41; p. 13 



Guo et al. 



THEMIS/MTR Fe 6302.5 A Bz 23-May-2008 17:12:29.450 UT 




500 510 
X (arcsec) 



THEMIS/MTR Na 5896±0.08 A Bz 23-May-2008 17:12:29.450 UT 





X (orcsec) 

SOT/SP Fe 6302,5 A Bz 23-May-2008 17:12:29,450 [U 




480 49 



X (arcsec) 

THFMIS/MTR Ha 6563 A 23-May-2008 17:12:29,450 UT 




X (orcsec) 

Figure 2. LOS magnetic fields of THEMIS/MTR Fe I 6302.5 A, Na Di 5896±0.8 A, and 
SOT/SP Fe I 6302.5 A in the solid rectangle region as shown in Figure [T] The Ha intensity 
map is shown in the bottom panel. Some slices are drawn and denoted by the solid lines on 
the magnetic field images, which will be used in Figures [6] and [S] 
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Figure 3. The vector magnetic fields of THEMIS/MTR Ca I 6102.7 A, Fe I 6302.5 A, Fe I 
5250.2 A, and SOT/SP Fe I 6302.5 A in tlie solid rectangle region as sliown in Figure[T] The 
180° ambiguity has been removed by the non-potential magnetic field calculation method and 
the fields have been transformed into the heliographic coordinates. Backgrounds are the vertical 
fields, white (black) is positive (negative). Arrows arc the horizontal fields. The threshold of 
the plotted arrows is 50 gauss. The circles in the bottom panel denote the regions of divergence 
and convergence of the arrows. 
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Figure 4. Histograms of the product of field strengths and filling factors, the azimuth, and 
inclination angles in the regions as shown in Figure [2] in the heliographic coordinates system. 
Only the fields where the field strength is greater than 50 gauss are included. The azimuth 
angles are measured counterclockwise from the north direction viewed towards the Sun and the 
180° ambiguity has been removed. The inclination angles are referred to the vertical direction 
of the solar surface. 0° is away from the Sun center and 180° is towards the Sun center. 
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Figure 5. Tlie scatter plots show tlie comparisons of S^, Bx, and By obtained by UNNOFIT 
and MELANIE inversions and SOT/SP observation in the field of view of the solid rectangle 
shown in Figure [T] 
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Figure 6. LOS magnetic field distribution along the white cut (left column) and along the 
black cut (right column) observed by SOT/SP (solid line) and MTR (dashed line). Top panels 
display the fields obtained by the inversion of the original SOT/SP observed profiles, while in 
the bottom row the Stokes profiles were convolved by a 3" X 3" 2D boxcar function at each 
wavelength before inversion. See text for details. 
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Fe 6302.5 A vs. Fe 5250.2 A B^^s Fe 6302.5 A vs. Ca 6102.7 A B, 
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Figure 7. Scatter plots of the LOS magnetic fields obtained by THEMIS/MTR Fe 6302.5 A 
vs. Fe 5250.2 A (left panel) and by MTR Fe 6302.5 A vs. Ca 6102.7 A (right panel). The data 
points are in the field of view of the solid rectangle shown in Figure fTI 



Table 2. The free parameters fitted in UNNOFIT and MELANIE. 



UNNOFIT 


MELANIE 


Zeeman splitting AAh 


Field strength B 


Inclination angle tp 


Inclination angle i/> 


Azimuth angle if> 


Azimuth angle if> 


Line strength r]o 


Lino strength rjo 


Doppler width AAd 


Doppler width AAd 


Damping parameter 7 


Damping constant a 


Line wavelength Aq 


Doppler velocity v 


Parameter of source function b 


Source function gradient Bi 


Filling factor a 


Stray light fraction / 




Macro-turbulence velocity Vm 



Hasan, S.S., van Ballcgooijen, A. A.: 2008, Astrophys. J. 680, 1542. 
Holweger, H., Miiller, E.A.: 1974, Solar Phys. 39, 19. 



Table 3. Formation heights at the line center of Fe I 6302.5 A and Fe I 5250.2 A. 
The heights are measured from rsooo = 1. 



Method 


Formation height (km) 


Atmosphere model 




Fe I 6302.5 A Fc i 5250.2 A 




Bommier et al. (2007) 


262 210 


Maltby et al. (1986) 


Bruls et al. (1991) 


232 259 


VAL-C 


Sheminova (1998) 


404 440 


Holweger-Miiller 
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Figure 8. LOS magnetic fields along the white (left column) and black (right column) cuts 
as shown in Figure [2] Top row shows the distributions for Ca I 6102.7 A (dotted), Fe I 6302.5 
A (dashed), and Fe I 5250.2 A (dash-dotted). Bottom row shows for Na Di 5896 A with AA = 
0.08 A (solid), 0.16 A (dotted), 0.24 A (dashed), and 0.32 A (dash-dotted). 
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Figure 9. The sketch of a flux tube model and the parameters for the flux tube, /m, Vm, and 
Bm are the Stokes parameters and LOS magnetic field at the flux tube center. /, V, and B 
are the observed parameters convolved with the quiet sun region. The intensity of the quiet 
sun is /q. The half-widths of the tube and the spatial resolution are denoted by w^x and s, 
respectively. All the parameters except s depend on the wavelength position AA from the line 
center of Na Di 5896 A. Left and right panels show observations at the edge and at the center 
of the flux tubes, respectively. 
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